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CONSTRAINTS ON VELOCITY ANISOTROPY OF SPHERICAL SYSTEMS 
WITH SEPARABLE AUGMENTED DENSITIES 

J. An 1 
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ABSTRACT 

If the augmented density of a spherical anisotropic system is assumed to be multiplicatively separable to func- 
tions of the potential and the radius, the radial function, which can be completely specified by the behavior of 
the anisotropy parameter alone, also fixes the anisotropic ratios of every higher-order velocity moment. It is 
inferred from this that the non-negativity of the distribution function necessarily limits the allowed behaviors 
of the radial function. This restriction is translated into the constraints on the behavior of the anisotropy pa- 
rameter. We find that not all radial variations of the anisotropy parameter satisfy these constraints and thus that 
there exist anisotropy profiles that cannot be consistent with any separable augmented density. 

Subject headings: galaxies: kinematics and dynamics — methods: analytical — dark matter 
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1. INTRODUCTION 

One basic problem in stellar dynamics is to find the dis- 
tribution function that is consistent with the given local den- 
sity profile. For the simplest cases, the problem reduces to 
solving an integral equation. The solution to finding the er- 
godic distributio n function f or an isotropic spherical system, 
i.e., the so-called Eddington formula is known to as early as its 
namesake. On the other hand, the complete solutions to devis- 
ing the two-integral even distribution function for an a xisym- 
metric system are also avai lable through the works by F ricke 
dl952l) . lLvnden-Belll ( ll962l) . and lHunter & Oianl dl993h etc. 

By contrast, the construction of a two-integral distribution 
function for an anisotropic spherical system contains an addi- 
tional difficulty. This is because integrating the two-integral 
distribution function over the velocity space produces an aug- 
mented density, which is a bivariate function of the potential 
and the radius. While it is an easy exercise to demonstrate 
that deducing the distribution function from the augmented 
density is a formally identical problem to the case of the ax- 
isymmetric s ystem with a two-inte gral even distribution func- 
tion (see e.g.. lOianll993tlAnl201 ll) . it is also obvious from the 
onset that, given the local density and the potential, there is no 
unique way to specify the bivariate augmented density with- 
out any additional prescription regarding the system. In fact, 
if the potential is known, specifying the augmented density 
is essentially equivalent to know ing an infinite subset of th e 
velocity moment functions (cf., [Dejonghe & Merrittl 119921) . 
to the zeroth of which the local density corresponds. The 
usual line of attack is to restrict either the augmented density 
or the distri bution function to be in the specified functional 
form (e.g.. lOsipkovl [1971 iMerritJ fl98l iCuddefordl U99ll 
Cudde ford & Louislll995t lAn & Evans! I2006at IWojtak et al. 
20081) and match further properties such as the anisotropy pa- 
rameter to that following the particular ansatz. One advan- 
tage of this approach is that these procedures usually simplify 
the subsequent inversion for the distribution function from the 
augmented density although these tend to sacrifice the flexi- 
bility in the behaviors of the varying anisotropy. 

On the other hand, the procedure that uniquely speci- 
fies the system and also allows the greatest possible free- 
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do m for the radia ll y vary ing anisotropy has been outlined 
in lOian & Hunted (11995b . If the potential- and radial- 
dependences of the augmented density are assumed to be mul- 
tiplicatively separable, then the radial part can be determined 
by the anisotropy parameter alone. The potential part then fol- 
lows immediately once the potential and the density are spec- 
ified. In fact, demanding the separable augmented density is 
the only route that fixes the anisotropy parameter, indepen- 
dently of the knowledge on the potential and the density. A 
practical implementation utilizing a general parametric form 
of the monotonica l ly vary ing anisotropy parameter is found in 
iBaes & Van Hesd ((2007), although adopting any parametric 
form intrinsically restricts the full flexibility of the technique. 

Recently, the separable augmented density has been at- 
tracted renewed interests in the context of the necessary con- 
ditions f or the distribution function to be non-negative. For 
instance, Ciotti & Morganti (2010b) proposed, for any system 
with a separable augmented density, the existence of the so- 
called global density slope-anisotropy inequality, which con- 
dition ally extends the centr al density slope-anisotropy theo- 
rem of lAn & Evans! d2006bl) to all radii. This has been subse- 
quently proved with a rest riction on t he central anisotropy by 
IVan Hese et all (l20Tlh and lAnl d201lh . 

This paper further explores the implications of the separa- 
bility in the augmented density and its limitation. In partic- 
ular, we examine the relation between the distribution func- 
tion and the moment functions and derive how the anisotropy 
parameter together with the potential and density uniquely 
specifies the separable augmented density. We also find that 
the radial part of the separable augmented density completely 
determines the anisotropic behaviors of not only the velocity 
dispersions but also every higher-order even velocity moment. 
Based on these findings, we also show that the non-negativity 
of the moment functions, which follows the non-negativity of 
the distribution function, restricts the physically permitted be- 
haviors of the separable augmented density and the anisotropy 
parameter described by it. In the following discussion, we 
consider our constraints in the context of the necessary and 
sufficient conditions for a separable augmented density corre- 
sponding to a physical system, extending prior works. Finally, 
further discussion concerning the inversion for the anisotropy 
profile in relation to the separable augmented density and our 
constraints is also provided. 
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2. PRELIMINARY 



The Ueansl theorem indicates that a steady-state spherical 
dynamical system is described by the phase-space distribution 
function (DF) of the form of /(£, If), Here £, the specific 
binding energy, and L, the magnitude of the specific angular 
momentum, are the two classical isotropic isolating integrals 
admitted by the spherical potential, i.e., <5 = *P(r) - jv 2 and 

L = (L-L)2 = ru t where is the relative potential with respect 
to the boundary (Hence, & > for a bound particle). Finally, 
v = (v 2 + v 2 )'- with u t = (v 2 + v 2 )? and v r being the tangential 
and radial velocities. Note that (v r , vg, v$) constitutes the set 
of three orthogonal velocity components in a common unit. 
Integrating the DF over the velocity space, 



NOP, r 2 ) 



.01 
= 5Eff 



d 3 vf(6, L 2 ) 



fd&dL 2 



(1) 
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results in a bivariate function of and r 2 . This is usually 
referred to as the "augmented density" (AD). Once the po- 
tential *P = *P(r) (which is not necessarily generated by the 
following density) is specified, the local density is found to 
be v(r) = N[*P(r), r 2 ]. In a self-consistent system on the other 
hand, the Poisson equation with the AD as the source term 
results in an ordinary differential equation on T(r), which can 
be solved to determine v(r) uniquely. 

The local higher-order velocity moments are also found 
similarly. Whereas any odd-integral moment must vanish 
thanks to the spherical symmetry (in particular, the isotropy 
in the configuration space), all the even-integral moments are 

found to be vf'v^ = m Pi? /N where 

JJJv 2 <XV 



In 



jj d&dL 2 <K p ^-L 2q f. 



(2a) 



Here the transform kernel is given by 

L 2 

9C(S, L 2 ; r 2 ) = 2(¥ - S) - - T , 



(2b) 



which is actually u 2 expressed as a function of the 4-tuple 



(fi, L 2 ; y ¥, r 2 ), whereas the integral is over the region in (£, L z ) 
space defined to be 

rs{(£,L 2 )|£>0,L 2 >0,7C>0), (2c) 

that is, the triangular region bounded by lines fi = 0, L 2 — 
and 7C = 0. The last line is the same as the diagonal line given 
byfi + (2r 2 )-'L 2 = ¥. 
If p > i, we find that 



<9(r 2 « +2 m M ) 



<9X 



(2p - 1) Trjjrffi c/L 2 *H ^L 2 "f (3) 



where X = or r z . Given that = 2 and 



W ~~ " SO 2 ) ~ 7^ eC l Ua " 
tion (0) indicates the existence of differential recursion rela- 



tions for the moment functions 



dm,, 



<9¥ 

d(r 2 * +2 m M ) 
dr 2 



(2p — l)m 



p-l,q, 



= {p-\)r 2 im p - 



1,4+1 > 



(4a) 
(4b) 



which is valid for p > i. In fact, once the AD is specified, 
every other velocity moment can be recovered without invert- 
ing for the DF. Specifically, given mo,o = N and the 'initial 
conditions' o(0, r 2 ) - 0, we first find that, for k > 1 



k times 



m h0 Q¥, r 2 ) = 2 k {\\ dW k - ■ • tWi moflQVu 
Jo Jo 
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r 2 ) 



(5a) 



g(g)j 



via repeated integrations of equation (|4al> and the Cauchy 
formula for repeated integration (eq. IA5ct . Here (a) n = 
n"=i( a + i~ 1) is the rising sequential product (the Pochham- 
mer symbol). Next, the repetitions of equation (l4bt yield 



i> 



-j) 



- ^ ('"Vf ( '" 2mt ' 0) ' (5b) 



for < q < k. Com bining these and using eq uation dA3K we 
recover the result of Dejonghe & Merritt (1992, eq. 13), 



m k - m Q¥,r 2 ) 



(k-iy. 



j"dQ Qf - Qf-^Ji^NiQ, r 2 )] (5c) 



for k > 1 and < q < k. The proper verification that 
equation d5cl ) is the unique solution to equations © given 
m p> q(0, r 2 ) — is provided in Appendix [A] An immediate 
corollary following equation (l5ct is that 



(A)\^N(T, r 2 )]>0 



(6) 



for every non-negative integer q is a sufficient (but not neces- 
sary) condition for every m PA to be non-negative. 

The behavior of the anisotropic velocity dispersions in a 
spherical system is u sually param etrized by the "velocity 
anisotropy parameter ' (Binn ey|[T980h . 



or 



2m lfi mr),r 2 ] 



(7a) 



where cr 2 = v 2 and cr\ = v\ = u 2 + u 2 . Meanwhile, equa- 
tion d4bl with (p, q) = (1,0) reduces to 



d{r 2 m\ fi ) 
dr 2 



2 ' 



(7b) 



Hence, the anisotropy parameter is directly related to the ra- 
dial partial derivative of the moment function m\#, i.e., 



0=1- 



2m, o 



5 In mi o 



d\nr 2 



V(r),r 2 



(7c) 
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Then, the total radial derivative of toi o results in 
dm u) 



dr dr + d*P ~dr~ 

2mio51nmio d^V 

■m,o- 



(8) 



r Slnr 2 dr 

With *P = *P(r) and equation ( TTcl i. equation © is simply the 
second-order steady-state spherical Jeans equation. Note in 
fact that *P(r) = Oo - O(r) is the relative potential with re- 
spect to the reference <t>o where €>(r) is the true gravitational 
potential, and thus — ^ = ^ = -V 1 where M r is the enclosed 
gravitating mass within the sphere of radius r. 



More generally, for p > ~, equation d4bb indicates that 



dm D 



[(q + l)m M - (p - 



dr r 

The corresponding total radial derivatives result in 



(9) 



dm„ 



dr 



— [(q + l)m p>q - (p - i)m p _i i9+ i] 



d¥ 



+ (2p-\)m p - Ul —, (10) 

which in fact constitute the co mplete set of the Jeans equa- 
tions ( Deionghe & Merrittl 19921) - see also lMerrifield & Kenl 
(1990) for the fourth-order equations, which correspond to 
(p,q) = (2,0) and (1,1) here. 

3. SEPARABLE AUGMENTED DENSITY 

Let us suppose that the *P- and r 2 -dependences of the AD 
are multiplicatively separable as in 



NOP, r 2 ) = P(^)R(r 2 ), 



(11) 



for some PQ¥) and R(r 2 ). It then follows from equation ( f5cb 
that every moment function m p$g Q¥, r 2 ) is also separable. In 
particular, m p , q (W, r 2 ) = 2P + H\) p P p+q Q¥)R q (r 2 ) or 



m k ^„Q¥,r 2 ) = 2 k ( l 1 ) k _ n P k ^>)R n (r 2 ) 



for < n < k, where 

(PQf) 
W) = I l 



(k = 0) 



dQ^-Qf-'PiQ) (*>!) 



(12a) 



(12b) 



c) 



R " {rl) S (d^J [r2 " Rir2)] = ^ ( r4 d^J [r2R(r2)l (12 

We find from equation ( 112bl > that -jkP_k _ Pk-i f° r anv P os ~ 
itive integer k. Similarly, equations d!2ct and (lA3t lead to 
^j(r 2n+2 i?n) = r 2 "R n +\ for any non-negative integer n. Next, 
since mioOP, J" 2 ) = ^ > i(*P)^('" 2 ), equation ( T7cT > indicates 
d In R(r 2 ) 



Pir) = 

R(r 2 ) 
R(f 2 ) 



d\n r 2 
exp 



(13a) 
(13b) 



In other words, the radial function R(r 2 ) is completely speci- 
fied (up to an immaterial scale constant) given the anisotropy 
parameter f3(r). Once R is specified, the potential part imme- 
diately follows the local density as PQ¥) = v(r )/R(r 2 ) with the 
inverse function r — f" 1 (*P) o f the potential (Oian & Hunted 
[19951: iBaes & Van Hesdl2007h . 



3.1. Implications of the separable augmented density 

Given the boundary conditions *P(ro) = and wi o(0, r?) = 
at r — ro (which may be the infinity), the radial velocity 
dispersion is given by 



dQP(Q) (14 

o 



R(r 2 ) 



v(r) d*¥ 



R(r 2 ) dr 



a) 
(14b) 



Here equation (I14bb is actually the solution to the stead y-state 
spherical Jeans equation with j?" 1 of equation d 1 3bt being 
its in tegrating factor (e.g., Ivan der Marellll99"4l lAn & Evansl 
2009). In other words, equation (I14bt always provides the ve- 
locity dispersions of the system given the potential, the den- 
sity, and the anisotropy parameter irrespective of the separa- 
bility assumption. 

However, the true implications of the separability assump- 
tion on the other hand lie beyond the behaviors of the velocity 
dispersions. That is to say, with the separable AD assump- 
tion, the anisotropy parameter not only specifies the complete 
AD together with the local density and the potential, but also 
it constrains the anisotropic behaviors of every higher-order 
velocity moment (including naturally those of velocity dis- 
persions) by itself. In particular, 



_ "« ,\ , \ m P,n 

°« = F = ( 2 + m ' 

«0 Wlp+nX) 



(15a) 



while a n (r) for any non-negative integer n is determined re- 
cursively from y6(r) alone (AppendixlBl via 



Jfn+i =(n + l -P)a n + 



with at) — 1 . Here, a' = = 
this is specifically translated into 



r da n 
2~dr 



a' n (15b) 
For a few small n's, 



a 2 



(1 



3 m 



1.1 



5m 2 - 



2m i,o 



2m 2 fi 2to 3 ,o 

3mo2 15mi 2 



4m2,o 4m3.o 
a 3 = (1 -j8)(2-j8)(3 -0 - 3(2 -/?)/?' -j8" 



15mp.3 
8m3,o 



(15c) 
(15d) 

(15e) 



and so on. Here, /?" = = t~ where u 
Furthermore, we also have 



lnr 2 etc. 



dm„ 



dr 



din 
dr\m 



p+qfi 
C Pl qRq 



1p+q,0J 



dm 



p+q,0 



c 



dm„ 



dr 



dr 
d tR q 

- C P R l d Il + c 



+ Cp,qPp+qPo~ 



da q 
dr 



p,qPp+qPo~ 



da q 
dr 



dm 



p,q 



5*P 



= C 



dP 



p,g 



P +C I p _ r p n 

^vp 1 ~ ^p,q r p+q-^q 



(16a) 

(16b) 
(16c) 



where C PA = 2 p+q ( j) p . Therefore, expressing the total radial 
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derivative of m Piq QP, r 2 ) leads to (note Rq = R) 



dm 



P+1,0 



dr 



= C 



p+qfi 



p+q 



dR 
~dr~ 



+ P 



p+q-l 



~d7 



R- 



(17a) 



Given equation ( 113al ). this indicates that, if the AD is assumed 
to be separable, the set of the (2n)-th order spherical Jeans 
equations (eg. [TOl with p + q - n and p > 1) reduces to a 
single equation, 



dm nQ 2f3 d*¥ 
— — ^- + — m„ ;0 = (2n - l)m„_i,o— . 
dr r dr 



(17b) 



This generalizes the fourth-order Jeans e q uation for 'con- 
stant anisotropy' intr oduced by lLokasI (120021 see also 
lLokas &"M amon 2003). We note however that equation ( I17bb 
is in fact the result of the separability assumption and not of 
the constant anisotropy per se. In addition, under the sepa- 
rability assumption, the solution to equation ( 117bt is imme- 
diately obvious as per m„o = (2n - l)!!P„[ v F(r)]/?(r 2 ) with 
equations <fl2bl) and ( fT3bl as well as P\^(f)\ = v(r)/R(r 2 ). 

3.2. Constraints on the anisotropy parameter 
3.2. 1 . general cases with a separable augmented density 

The non-negativity of the DF implies that all the even- 
integral moment functions must be also non-negative. Conse- 
quently, if N(T\ r 1 ) = PC¥)R(r 2 ) is separable, then P n (W) > 
and R,,(r 2 ) > for any non-negative integer n. While PQ¥) > 
is the sufficient (and also the necessary since P = Pq) con- 
dition for P„(*F) > for any non-negative integer n, the con- 
dition that 



, d"\x n R(x)] 



dx n 

or equivalently (see eq. |A31 > 



> (r, x > 0) 



(18a) 



(-1)" 



d" Riw- 1 ) 



dw" 



> (r, w > 0) 
(18b) 



for every non-negative integer n constitutes a set of indepen- 
dent constraints on the behavior of /?(r 2 )Q In other words, 
equation ( fT8l is a necessary condition for the radial part R(r 2 ) 
of any separable AD to be g enera ted by a non-negative DF. 

Combined with equation d!3at , it forms the set of restric- 
tions on the radial variations of f}(r) allowed for the spheri- 
cal system with a separable AD. In particular, any spherical 
anisotropic system with a separable AD is physical only if 



a n (r)>0 (r>0) 



(19a) 



for all positive integ ers n . Here the set of functions a„ is as 
defined in equation ( I15bl ) with f3(r). For the first few small 

' A function <p(x) of x > is said to be "completely monotonic" (cm.) if 
and only if (-l)"0 ( "'(x) > for all non-negative integers n. Hence, the condi- 
tion is equivalent to saying K(w) = R(uT )/wis a c m. function of w. Accord- 
ing to S. Bernstein's theorem on cm. functions (see Widder 1941), an impor- 
tant corollary to this is that K(w) must be the Laplace transformation of a non- 
negative function - i.e., the inverse Laplace transformation of K(w) exists and 
is non-negative for positive reals. Moreover, expressing the inverse Laplace 
transf ormation using E. Pos t's inversion formula (see Hirschman & Widdeil 
1955) indicates that eq. j I S I is actually equivalent to lim,,-,,*, R n ( j)/«! > 
for t > 0. Further explorations of this idea will be given elsewhere. 



n's, the conditions are equivalent to 

p< i; 



r dB 

-f r <{\~P){2-P)\ 



(19b) 
(19c) 



r d , dB\ 3 dB 



and so on. Equation ( 119b| ) is obvious from the definition of 
the anisotropy parameter (eq. |7a) and the non-negativity of the 
velocity dispersions, and thus universal independently of the 
separability assumption. By contrast, the further constraints 
involving the radial derivatives of ft are the consequence of 
the separability assumption - following equation (TS[ and the 
non-negativity of the higher-order velocity moments. 

These imply that, even if R(r 2 ) could be formally written 
down using equation d!3b| ), not all arbitrarily varying f3(r) are 
consistent with separable AD because some might produce 
negative higher-order moments. For example, consider the 
anisotropy parameter behaving 



„2.v 



/J(r) = 



(20a) 



so that the system is isotropic at the center and radially 
anisotropic in the outskirts. I f s — 1, the an isotropy profile 
of equation d20a| ) is that of the Osipkov- Merritt (OM) system. 
The corresponding radial function is R(r 2 ) = (1 + r 2s jr 2s )' l l s 
within a constant, but this does not satisfy the condition in 
equation ( fT8b if s > 1 because 



d 2 [x 2 R(x)] 



dx 2 



x=r 2 /rl 



2-(s-l)x s 
(1 +x s )^ +2 



(20b) 



which is negative for x > ^J2/(s - 1). Equivalently we find 

2r 2s -(s- \)r 2s 



2s\2 



(20c) 



(rf + H-0 

and thus equation J20ab fails the constraint in equation (I19cb 
if s > 1 a nd rj r a > H]2l{s - 1). Consequently, /3(r) in equa- 
tion (120al i is consistent with a separable AD only if s < 1. 
In fact, the converse is also true, i.e., if s < 1, then R(r 2 ) = 
(1 + r 2s lr 2s y l l s satisfies equation (TT8l for all non-negative in- 
tegers n - obviously, if s — 1, the non-negative OM DF exists 
with a properly chose n pot ential term ^(^10. 

Roughly, equation d!9cb insists that the anisotropy param- 
eter in the system with a separable AD cannot increase ra- 
dially faster than the limiting value determined by the local 
anisotropy parameter, which tends to get smaller as it be- 
comes more radially anisotropic. Similar interpretations for 
higher-order constraints of equation (|T9l are less obvious. 

3.2.2. a family of monotonic anisotropy parameters 
Consider the anisotropy profile, 



j8(r) = 



„2.s 



,2 s 



(s > 0). 



(21a) 



This parametriz a tion has also been introduced by 
iBaes & Van Hesd d2007f) for their construction of dynamical 
models with a flexible anisotropy parameter. Equation (120a) 
corresponds to equation (12 1 ab with (J3q,/3co) = (0, 1). Note 
that the transform of s — > -s is actually equivalent to 
switching /3o <-> and thus the restriction s > is 
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actually not necessary. Nevertheless, to assign definite 
physical meanings to the parameters, we retain the restric- 
tion. Then B monotonically varies from Bq at the center 
to /Joe as r — > oo, with the constant-/? case represented by 
Bq = jSoo- The choice that r a = or r a = oo also produces 
the constant-/? model although they will not be considered 
explicitly here. T he case s = \ reduces t o the generalized 



fails if Bq < ftca and 



Mamo n & Lokasl anis otropy model (cf„ IMamon & Boud 
20101) with the original IMamon & Lokasl (120051) model given 
by (/3o,0oc) = (0, j). The sepa rable AD with s = 1 incl ude 
the Cuddeford system (see also Ciotti & Morganti 2010a) for 
which Poo = 1 and the O M system with (fio,Boo) = (0, 1). As 
iBaes & Van Hesd (120071) have noticed, this parametrization is 
notable as it yields the simple anal ytic in tegrating factor for 
the second-order Jeans equation (eq. I13bh 



[R(r 2 )T l = r 2 ^(r 2s + r 2 d s ) 



(21b) 



With v(r) and *F(r) specified, the radial velo city dispersion 
can be found in quadrature by equation dl4b) . regardless of 
the separability of the AD. 

Under the separable AD assumption on the other hand, 
this completely specifies the resulting system; N( V F, r 2 ) = 
PQP)R(r) where P[¥(r)] = v(r)/R(r 2 ) and R{r 2 ) is de- 
duced from equation (I21bt . The DF can be found by in- 
vert ing the integral equatio n ([T) - e.g., Deionghe (11986b 
and IBaes & Van Hesd (l2007h for the tec hniq ue based on the 
Laplac e-Mellin transform or lOianl (11993b and iHunter & Oianl 
(I1993L see also lAnl 1201 ll) for the complex contour integral 
method. 

However, the preceding arguments indicate that the result- 
ing model is not necessarily physical for an arbitrary param- 
eter set. Obviously, the condition that B{r) < 1 for Vr > 
restricts the parameters to be Bq < 1 and Boo < 1. As the 
system with a separable AD, more constraints on the param- 
eters also follow equations (fT8l and ( fT9l . The first of these 
corresponding to equation (fT8l with n — 2 reduces to 

[(2 -j8 ) + (2 -Poo)x s ] [(1 -Bq) + (1 -Poo)x s ] 

>s(fioo- Bq) x s (for x > 0), (22a) 

which is also equivalent to equation ( |19ct ; 



(r 2s + ly 



= [(2 -/3oo)r 2s + (2 -M [(1 -Poo)r 2s + (1 - Bq)} 

-s(0oo- Bq) r 2s > (for r > 0). (22b) 

Here, we have set r. d — 1 for brevity, but this does not affect the 
following results. Since the necessary and sufficient condition 
for the real-coefficient monic quadratic equation x 2 +bx+c = 
to possess no non-degenerate positive real root is c > and 
b > -2|c|i, equation d2"2l > for Bq, Boo < 1 is also equivalent to 



(2 - A,)(l - /Jo) + (2 - B )(l - Poo) - s(J3oo - A>) 

+ 2(2 - j8b) (2-^)1(1-^)' >0. (22c) 

With(2-/U(l-A)+(2-A))(l-£co) = 2(2-A))(l-y6co)+ySoo- 
jSo> we therefore find for fixed fio,fioo < 1 that equation d2"2i i 



s > 1 + 



2(2-^,)s(l -pooY 



x [(2-#>)*(l -^r- + (1 -Po)H2 -#»)*]• (23) 

That is, there exist parameter combinations for equation (|2T1 > 
that cannot be consistent with any physical separable AD. 

More higher-order constraints may be derived similarly, but 
direct calculations for general cases become rather compli- 
cated as the order increases. Instead, here we just note that, 
if < s < 1 o r fo'P™ is a non-negative integer, then Rir 2 ) in 
equation d21b| i satisfies the condition of equation ( fT8l and so 
yS(r) in equation (121 at can be consistent with a physical sepa- 
rable AD. An elementary proof is provided in Appendix |C] 
The sufficiency of the condition that < s < 1 for the 
parametrization given in equation (|2H can also be deduced 
by the existence of the correspond ing non-negative DF wit h 
a separable AD as demonstrated by Bae s & Van Hesd (120071) . 
who explicitly constructed the particular DF in terms of the 
convergent Fox H-function. We also suspect that if fio < A», 
the condi tion t hat < s < 1 is the necessary condition for 
equation ( 12 lbb to satisfy equation ([T"8l but have no definite 
proof at this time. 

In addition, we also note that the condition in equation ( TT8l 
is linear on R{r 2 ). Hence, if both A(r 2 ) and B(r 2 ) meet the nec- 
essary condition in equation (fT8l l. the radial function given by 
the linear combination, R(r 2 ) = aA(r 2 ) + bB(r 2 ) where a and 
b are positive constants, also satisfies the same necessary con- 
dition (and therefore the anisotropy parameter resulting from 
it is consistent with eq.[l9]and a physical separable AD). For 
example, th is indicates that th e multicomponent generali zed 
Cuddeford systems studied by ICiotti"& Morganti (201 Pah do 
satisfy equation ( fT8l as their radial functions are g iven by the 
sums of the functions in the form of equation (12 lbb with 5=1 
(and Boo = 1) and different r a 's. 

4. DISCUSSION 

4. 1 . Sufficient conditions for physical separable augmented 

densities 

For multicomponent Cuddeford systems ICiotti & Mo rganti 
d2010al) have proved that the condition > is sufficient 
to guarantee the non-negativity of the posited DF. Here p, is 
the integer floor of (i.e., the greatest integer not larger than) 
| —Bq. Subsequently. IVan Hese et ail d201 lb asked whether 
the same condition should be the sufficient condition for any 
separable AD to be generated by the non-negative DF. 

The present paper clarifies the answer to their question in 
the simplest form to be negative. The existence of the nec- 
essary condition involving only the radial function R(r 2 ), i.e., 
equation ( fT8l . implies that any sufficient condition must also 
contain some restrictions on the same. The hypothesis as 
stated entirely with the potential part PQV) thus cannot be a 
sufficient condition by itself given the independent nature of 
the potential and radial parts of the separable AD. A sim ple 
counterexample may be constructed with equation (12 lab and 
the choice of parameters such that | < Bq < Boo — 1 and 
s > 1. With the radial function R(r 2 ) that follows (eq. l21bb . no 
function P(*P), regardless of ^ > or not, can lead to a non- 
negative DF - here the choice of P( V F) is equivalent to speci- 
fying the local density as P[*¥(r)] = v(r)/R(r 2 ) and insisting 
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^ > im poses the so-called global density slope anisotrop y 
inequality dCiotti & Morgantill20105 IVanHese et"aT1l201 ll) . 
Nonetheless, the findings in this paper do not preclude the 

possibility that the constraint ^ combined with ad- 

ditional conditions on R(r 2 ) may constitute a sufficient con- 
dition for the system with a separable AD. In fact, after the 
original version of this paper was submitted, E. Van Hese (pri- 
vate communication) has discovered the existence of the set 
of such conditions, e.g., together ^ > and the set of condi- 
tions on R(r 2 ) that includes equation ([T8l > are sufficient for the 
existence of a non-negative DF. 

4.2. Universal constraints on the anisotropy parameter? 

Since the constraints in equation ([T9l are actually put on 
the anisotropy parameter without any explicit reference to the 
separable AD, it seems fair to ponder how general these con- 
straints actually are - i.e., whether the constraints exist for 
any physical AD. Although we find, for the time being, no 
reason to argue that these constraints are universal (for they 
are derived based on the particular assumption of the separa- 
ble AD) beyond the obvious restriction that fi < 1, we will 
not attempt to settle the answer in this paper. However, we do 
note that, if the constraints are entirely the consequence of the 
separability of the AD, one must be able to construct a pair of 
the non-negative DF and the inseparable AD for a spherical 
anisotropic system with its anisotropy parameter violating the 
conditions in equation ( [19) . Unfortunately, the task is com- 
plicated by the fact that, with inseparable AD, the anisotropy 
parameter cannot be specified independently without impos- 
ing the fixed behavior of the potential; the pair of /(£, L 2 ) 
and N(*P, r 2 ) typically prescribes only y6( v F, r 2 ) and thus vary- 
ing *P = ^(r) results in a different /3(r) = ^[^(r), r 2 ]unless 
¥(/-) and v(r) = Nf^r), r 2 ] are related to each other through 
the Poisson equation and so the freedom to choose *F(r) is 
subsequently removed. 

4.3. Separable augmented densities and the Jeans 
degeneracy 

With real data, our observations are typically limited by 
projection, and thus the usual kinematical observables avail- 
able to us are restricted to the surface brightness profile (or the 
column density profile for the discrete number count data) and 
the line-of-sight (los) velocity dispersion. It is a well-known 
fact that, while the three-dimensional density profile can be 
uniquely inverted from the surface density (the Abel transfor- 
mation) under the spherical symmetry assumption, the radial 
and tangential velocity dispersions, cr 2 (r) and cr 2 (r) - which 
are needed to find the potential through the Jeans equation - 
cannot be determined from the los velocity dispersion alone as 
they are degenerate in reproduc ing the observations of the last 
(Deio nghdll987l : lMerrittlll987l) unless the system is known 
to have isotropic velocity distributions. One may lift this so- 
called Jeans degeneracy by imposing additional constraints on 
the system comi ng from observations o r a p riori assumptions . 

For example. lMamon & Bouel (120101) and lWolf et all (120101) 
have shown that it is in general possible to find cr 2 that 
is consistent with the observed los velocity dispersion pro- 
file and any arbitr arily specified anisotropy parameter f3(r). 
lEvans etafl (120091) on the other hand demonstrated that one 
can also go in the other way around, finding /3 that is con- 
sistent with the observed los velocity dispersion and the arbi- 
trary assumed form of cr 2 . Alternatively, under the assump- 
tion of the constant mass-to-light ratio, one can also find 



the unique solution to the coupled Jeans -Poisson equations 
from the los veloc i ty dispersion profile dBinnev & M amon 
[1981 iTonrvl [l98l iBicknell et all [19891 This is equiva- 
lent to specifying the potential first and inverting the los ve- 
locity dispersion to determine given the Jeans equatio n 
(ISolanes & Salvador-Soil 119901: iDejonghe & Merrittl 09921) . 
However, for the purpose of constraining the potential, these 
anisotropy inversion algorithms can only, at best, reject some 
choices for the gravitational potential as unphysical where fi 
reaches values above unity (formally this indicates negative 
velocity dispersion). 

The best observational constraint for lifting the Jeans de- 
generacy on the other hand woul d be some handle on the 
proper motions of tracers (e.g ., Leonard & Merrittl 119891: 
Ivan der Marel & Andersonl2010l) as they are the velocity pro- 
jections that are orthogonal to the los velocity. Alternatively, 
with a data set consisting of discrete tracers, the precisely 
measured differential distances (which ultimately yield the 
distances to the center of the system ) can also break the Jeans 
degeneracy (cf., Wat kins et al.ll2010T) . Unfortunately, with our 
current and near-future observational capabilities, their uses 
are mostly limited to very near-by objects. 

A popular idea for possible ob servational constraints is the 
use of the higher-order moments (Merri field & Kenfl [T990) or 
the distribution of the los velocities - note that specifying the 
distribution is equivalent to knowing the infinite set of the en- 
tire moments. Simila r to the velocity dispersion s (which are 
the second moments), Deiong he & Merrittl (119 92) have shown 
that, with the potential specified, the complete set of indepen- 
dent velocity moments can be solved from the observed los 
velocity moments up to the same order - in the infinite order, 
this implies that the DF is uniquely specified by the observed 
distribution of the los velocities provided that the potential is 
known a priori. However, it is easy to argue that this will not 
solve the degeneracy problem (in particular for tracing the po- 
tential observationally) because, under the spherical symme- 
try, introducing each new (2«)-th moment adds (n + 1) new 
variables and the n constraining Jeans equations (eq. IT0T > with 
one further observational constraint into the mix and therefore 
th ere is no net in creas e in the constraint s . 

iLokasi (120021) and lLokas & Mamonl (120031) on the other 
hand introduced a hybrid of theoretical and observational con- 
straints, 'constant anisotropy' and the fourth moments (kur- 
tosis), to bring the degeneracy problem into the unique so- 
lution. However, their 'constant anisotropy' is actually in 
the form of a strictly stronger assumption that the DF and 
the AD are given by the an satfl /(£, L 2 ) = L~ 2/1 g(S) and 
N(*P, r 2 ) = r~^ P(*F) while § [XT] of the present paper (in par- 
ticular, eq. H7bl > indicates that this is a rather unnecessarily re- 
strictive assumption for their method to work. That is to say, 
under the separability assumption of the AD, to bring higher- 
than-the-second-order moments into th e pro blem only adds 
one independent n ew va riable (cf., eq. 1 15 ab and the single 
Jeans equation (eq. I17bt . Therefore the observations of the 
los velocity moment at the same order actually act as an ad- 
ditional net constraint on the system given the separable AD. 

2 Recently, Woitak et al. 1 2008) proposed an extension of this by introduc- 
ing a more general ansatz for the L part of the DF that allows the variation of 
ft, which they found to be consis tent with the DF of simulated ACDM halos. 
In a sequel, Woitak et al. 1 2009) predicted the distribution of the los veloci- 
ties with which the observed data can be fit to determine the parameters of 
the DF. In principle, this is still less flexible than the procedure outlined in the 
following, but if in practice one were to parametrize the anisotropy parameter 
in a particular functional form, the approach may be seen as complementary. 
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Specifically, the introduction of the fourth moment is enough 
to uniquely solve the Jeans degeneracy if the AD is assumed 
to be separable whereas adding further higher-order moments 
actually over-constrains the problem. 

However, the assumption of a separable AD is purely for- 
mal and its physical interpretation is unclear, although it is a 
weaker hypothesis than the power-law ansatz for the L part 
of DF and the AD which produces the constant anisotropy. 
The most conservative statement that can be drawn regarding 
the Jeans degeneracy and the separable AD is thus that given 
the observations of the second and fourth moments (the dis- 
persion and the kurtosis) of the los velocities, there exists a 
unique spherical model with a separable AD that is consis- 
tent with them. The resulting model is complete in that it es- 
sentially specifies the DF as well as the underlying potential. 
Although the non-negativity of the resulting DF is not guar- 
anteed, the condition in equation ( TT~8T > is both necessary and 
sufficient to prove that the model will produce non-negative 
velocity moments of every order - of course, the model is not 
necessarily 'real' and the predicted higher- than-fourth- order 
moments should be compared to the observations (if avail- 
able) in order for it to be acceptable. 
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APPENDIX 
A. PROOF OF EQUATION gcjl 

Lemma Al For a non-negative integer n and arbitrary real a, 



d n (x*) 
dx" 



j=0 



(Al) 



which is easily proved by the induction on n. For a positive integer 
power monomial (i.e., £ is a non-negative integer), this simplifies 



| 77—77 (0<n<k) 



dx" 



(k-n)\ 



(Ala) 



(n>k+l) 
Next, using Lemma fATl and the extended Leibniz rule, we find that 
Lemma A2 for a non-negative integer n and any function f, 

r m-k{ ^\ jk r " 



d n {x»f) A jn\ d"- k (x") d k f 



z 



dx" £\k} dx-* „„ k __ Q 



*l-f-¥*—*f». (A2) 
dx k (k\y-(n-k)\ J 



Here M is the binomial coefficient and f (k \x) = Now, we are 
able to prove 



Theorem A3 for a non-negative integer n and any function f, 

(4)',/>-^- 



Proof. We prove this by the induction on n. First, equation dA3t is 
trivial for n = 0, 1. The induction step is proved as 



V ax I dx \ dx! 

(«!) 2 



dx 



dx" 



X dx[Z (*!,:,„ .. k)\ 



(k\) 2 (n-ky. 
(n\f[(n + 1 - k)(n + 1 + k) + k 2 ] „ 



(A3a) 



+* Ak) 



(k\) 2 (n-k+ 1)! 

= x «+2 y (n\) 2 (n+l) 2 

Z^( k \nn-k+l)\ ■> 



x" 



x n+K f 



+1 d" +l (x" +l f) 
dx" +l 



Immediately following this is 

Corollary A4 for any non-negative integer k, 

dxl dx k J dx[\ dx/ 

1 I d d k+l (x k+[ f) 



(A4) 



dx k+l 



With / = N( V P, r 2 ), x = r 2 , and k = q, this results in 

^[^(^ft^NOP.r 2 )]] = r 2 \^f\r 2 "^,r 2 )]. (A4a) 

It is now easy to show that m pr/ C¥, r 2 ) in equation (TScb satisfies 
equation d4bb by direct calculations using equation dA4at and the 
Pochhammer symbol (j) p = (^) p -i(p — |) for any positive integer p. 
One can show that equation i5c[ satisfies equation ( Hat using 

— f dij{x-y) k g(y) = S ka g(x) + k f dy (x - yf' 1 g(y) (A5a) 



(where d„,,„ is the Kronecker delta) and 



f dxi\ 

Jxn J.in 



dy (x\ - y) g(y) 



| dy g(y) I dx x (x x 



■yf 



(A5b) 



for any k > 0. By repeatedly applying this, one finds that a simple 
iterated integral in general reduces to an integral transform, 



9k(x) 



= f dx k - - f dx, g(x{) = 1 f dy(x- 
Jo Jo (fc-l)Uo 



dx" 



(A5c) 

where k is now a positive integer. This is sometimes known as 
the Cauchy formula for repeated integration, and can be strictly 
proved through the induction on k. In fact, the function g^(x) de- 
fined as such is the particular solution to the differential equation 
d k gi l (x)/dx k = g(x) with the set of initial co nditio ns ffjP(P) = 
where j e { 0, . . . , k — 1 ). Formally, equation &A5c\ extends to the 
case k = by noting that lim f ^ + e M e 1 = 2S(x) where S(x) is the 
Dirac delta. The extra factor of two is due to the fact that the integral 
interval in equation dA5ct is one-sided extending down from x. 
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B. GENERAL EXPRESSION FOR EQUATIONS O&O 
For any function f(x) of x, we find using equation JA4b that 



d" +[ [x" +l f(x)] 
dx" +1 



x" dx [ dx" \ 

[ l+n + x Tx) 



d \ d»[x»f(x)] 



(Bl) 



Next, with the definitions of a„ in equation J15at and of R„ in equa- 
tion l |12ct 

, . Kir 2 ) 1 d"[x-R(x)] 



R (r 2 ) R dx" 



(B2) 



from the given radial functi on R( r 2 ) of a separable AD, the recursion 
formula for a„ in equation d 1 5bl > follows as 



/ d\aR\ 

= ( l + n + 7m?r 



.v=r- 



(B3) 



dlnr 2 



Replacing R with ft utilizing equation j!3at yields equation dl5b|> . 
Equation J18t then indicates equation M9\ whereas equation d 1 5bl > 
may be considered to be the recursive definition of a n (r) from the 
anisotropy paramet er f}(r ) given the initial term a Q - i. Note also 
that a„ in equation J 15bt is defined as such without referring to the 
radial function or the separable AD at all. 

C. PROOF OF CONSISTENCY OF EQUATION {ZD 
For p(r) and R(x) given by equations l !21b . if we define 



(1 +x')"a„ =xP°(l+x s ) n+ 



dx" L(l + x s f 



(CI) 



where A = tea > equation d!5bt results in 



r n+ , = [(n + 1 -y8 ) + (fl + 1 -/?„ - sn)y]r„ + sy(l + y)-^, (C2) 

ay 

where y = x" = 5j. Since t = ao = 1, this indicates that r„ is an (at 
most) n-lh order polymonial of !/. If we then let 



t2T 



(C3) 



we can derive the recursion relation for the coefficients, 

?„+u = («+ 1 -/Jo + sfc)£a + tO - + 1 -#» + l)]?„,t-i (C4a) 

by substituting equation jC31 > into equat ion dC2t . and also using t n k = 
for k < or k > n. If k = 0, equation dC4at reduces to 



5w-i,o = (n+ 1 -fio%,c 



i n ,o = (1 > 



(C4b) 



provided that y6o < 1 because f„ _i = an d i = 1. Next, for a 
positive integer pair n > k > 1, equation dC4ab indicates that, if 
< j < 1, the non-negativity of f„_i_^ and f„_i_t_i can guarantee 
the non-negativity of t„^ provided that /?o,ySc» < 1. Since we have 
already found that = for k < and f„_o = (1 - y3o)„ > 0, we can 
conclude that if < s < 1 and f3o,/3 m < 1, then T„(y) is a polynomial 
with all non-negative coefficients and therefore T„(y) > for v i/ > 



and any non-negative integer n. Since t„ = (1 + x s )"a„ = (1 + x 5 )^, 
equation d 1 81 > also follows immediately. 

As for the cases that f = = -A is a non-negative integer, we 

first consider the constant-/? case, that is, £ = and /?o = = 
Then, using equation l lAlb . 



R n (x) ■■ 



d n x n-/3 

dx" 



for x > 0, provided that yS < 1. In general, if ^ = — is a non- 
negative integer, we can simply extend this result to 



(1 + x')t ri 



(C5b) 



'-2 



A d"(x sk+ "-k') ^ (g 



*=0 w i=0 



^mCsifc+l-A^Jt* 1 *. (C5c) 



Again, R„(x) > for x > 0, provided that j8«, = j6 - sg < fS < 1. 
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